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~i ■ The purpose of this paper is to calculate the two-photon decay rate corresponding 

to the two-photon transitions nS — > 15 and nD -^ \S in hydrogenlike ions with 
a low nuclear charge number Z (for principal quantum numbers n = 2, ... ,8). 
Numerical results are obtained within a nonrelativistic framework, and the results 

fT^ ■ are found to scale approximately as (Zoif' jv? , where a is the fine-structure 

constant. We also attempt to clarify a number of subtle issues regarding the 
treatment of the coherent, quasi-simultaneous emission of the two photons as 

^ rS ,' opposed to one-photon cascades. In particular, the gauge invariance of the decay 

rate is shown explicitly. 
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The subject of the current paper is the two-photon decay rate of excited atomic states, interpreted as the imaginary 
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O 1 Introduction 
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f^ ' part of the two-loop self-energy. We follow our previous investigation reported in Ref. [1] and augment the analysis 
by treating the decay rate in both length and velocity gauges. Special emphasis is placed on the role of singularities, 
infinitesimally displaced from the integration contours for the photon energy integrations, which are generated by 
' ;> ' bound-state poles of lower energy than the reference state (in the sense of the two-loop self-energy) . The reference 
QQ state is equivalent to the initial state of the two-photon decay process. A good quantitative understanding of the 
^^ ,: two-photon decay processes from highly excited hydrogenic bound states is important for astrophysics, as emphasized 
>0 ■ in a recent paper by Chluba and Sunyaev [2] . As the physics of the process is in principle well known and has been 
discussed in a previous fast- track communication [1], we see no obstacle to going in medias res with the analysis. 
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Our purpose here, in addition to providing numerical data concerning the D —> S transitions, is to clarify the role 
QQ ' of cascades of one-photon decays through so-called resonant intermediate states, which are addressed using concepts 
f^ . developed in field theory [3,4]. 

^ , Natural units with h — c — e^ — 1, i.e. e^ = Ana, are used throughout this paper, which is organized as follows. In 
Sec. [21 the gauge invariance of the two-photon decay rate, as derived from the two-loop self-energy, is reanalyzed. In 
Sec. [31 numerical results for nD — > IS' transitions are presented; these were not treated in the previous paper [1]. A 

d ' discussion of our results, including a comparison to previous investigations of two-photon decay from highly excited 
states (see Refs. [5-7]) is given in Sec. [4] Cascade contributions are analyzed in Sec. [5] Conclusions are drawn in Sec.[6l 

2 Gauge Invariance 

We start by considering the two-photon self-energy for a reference state \(j)i) in a hydrogenlike ion, as derived from 
nonrelativistic quantum electrodynamics (NRQED). In the velocity gauge, the interaction Hamiltonian of the quantized 
electromagnetic field with the electron is given by 

Hj^^—(p-A + A-p) + '-—, (1) 

2m V / 2m 

where A is the vector potential of the quantized electromagnetic field. 



The well-known expression (see, e.g., Refs. [1,8]) for the two- loop self-energy reads (wi and a;2 denote the energies of 
the two virtual quanta) 
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Here, ReAEl is the real part of the energy shift, which gives rise, in particular, to the so-called two-loop Bcthe 
logarithms [9]. Our treatment relies on the identification of the imaginary part of the energy shift in terms of a 
decay rate of the reference state, as suggested by Barbieri and Sucher in Ref. [10]. In Eq. ^, Sr\ ' is a correction 

(2) . 

to the one-photon decay rate, whereas F^- is the two-photon decay rate. The former is obtained by terms where the 

integration over uji or 0^2 meets a bound-state pole and generates an imaginary part, in the sense of Eq. (4) of Rcf. [1], 

(2) 
but the other photon energy is integrated with a principal- value prescription. The latter term, T^ , is obtained by 

selecting exclusively the imaginary part generated by the singularities at uji + uj2 = Ei — E^ , where E^ is a virtual 

state contained in one of the propagators. All expressions on the right-hand side of Eq. ^ are manifestly of order 

a^{Za)^m, i.e. {Za)^Roo where i?oo is the Rydberg constant. 

The function /^ reads as follows (with all infinitesimal imaginary parts duly taken into account), 
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where the terms denoted by the ellipsis are given in Eq. (3) of Ref. [1], being irrelevant for the current investigation, 
because the two-photon decay rate is generated exclusively by the poles where the sum uji + uj2 of both photon energies 
is on resonance. In a basis-set representation, the expression for the two-photon decay rate V^"^' is thus found from the 
first three terms in Eq. ^ as [1] 
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where we use the summation convention for the Cartesian coordinates labeled by the indices j G {1,2,3} and k € 
{1,2,3}. The sum over v contains all virtual states, i.e. over the entire bound and continuous spectrum. We here 
imply a sum over the magnetic projections of the intermediate states, and of the final state of the decay process, 
but an averaging over magnetic projections of the initial state (since the decay rate does not depend on the magnetic 
projection of the initial state, one may alternatively choose any allowed value for the initial-state magnetic projection). 

We now assume all initial and final, and virtual states to be given in terms of hydrogen wave functions in the standard 
representation (see, e.g., Ref. [11]), so that 

{(j)f Ip' I (j>y) ((/).„ Ip' I <j>i) = {(j)i Ip' I (/),„) {<j>y Ip' I (/)/) . (5) 

where the sum over j is assumed. We then do the angular algebra [12]. For nS —> 15 decays, one obtains a result [1] 
which reproduces the well-known expression obtained by Goppert-Mayer in Ref. [13] for the particular case of \cj)i) = 
|25), 

T^^^^^^irnRe^^ujujiE^S-Eis-uj) hl^''^^^^^'"''^ 
"^ 277rm2.^o J Vni 16 MZ^\ E„s - E,p - u, + ie Eis - E,p + uj + ie \ ' 

g \ I^ / 

(6) 
where we use the definition of the reduced matrix elements according to Ref. [12]. Virtual P states are also relevant 
for the decay nD -^ IS decays, but the well-known prefactor is different [7], and the result is 
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where for completeness we note that the reduced matrix element for P ^ D transitions differs from the "radial" 
component of the matrix element by a factor \/2- 

In the length gauge, the atom-field interaction is given by 

Hi = -eE-r, (8) 

where E is the quantized electric-field operator. The length-gauge two-photon self-energy is obtained by straightforward 
fourth-order perturbation theory as 
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We observe the factor wf wfi which is characteristic of the length-gauge formulation. The absence of the seagull term 
as opposed to the velocity gauge leads to a somewhat simplified expression, 
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In contrast to Eq. ([3]), the momentum operators are replaced by position operators. In analogy to Eq. ([3]), only the 
first three terms are relevant for the two-photon decay rate. Using a basis-set representation, the expression for the 
two-photon decay rate derived in the length gauge thus reads 
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where the sum over v contains all virtual states. Using the identity 
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it is easy to show the equivalence of the two expressions for the two-photon decay rate given in Eqs. ([2]) and ([9]). Note 
that this equivalence can be shown easily using the commutator relation p* — i [H, x^ ] , but it holds only if the sum 
over V extends over the complete spectrum. 

Assuming hydrogen wave functions in the standard representation, we have that in analogy to Eq. ([5]), 
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After angular algebra, one obtains for the decay nS -^ IS*, 
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whereas for nD ^ IS decays, 
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again in complete analogy to Eqs. JH) and ([7|), respectively. 
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3 Numerical Results 

We here focus on the nS -^ 15 and nD — > IS decays, as indicated in Eqs. (fT4|) and fTS]) . respectively. Decays to 
the ground state have the highest rate for both one-photon [14] as well as two-photon processes and are therefore of 
special interest. Due to the infinitesimal imaginary parts explicitly indicated in Eqs. (|14p and (|15p . we can extend the 
sum over intermediate, virtual states over the entire hydrogenic spectrum, including those P states which have a lower 
energy than the reference state. We recall here that the double poles at intermediate resonances are naturally treated 
using the formula [1] 

lim Re f duj ( 3__] ^^L_. (16) 
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Simple poles are treated using the well-known Dirac prescription, and the principal-value integration then yields the 
real part of the integrals. Numerical results can be obtained by expressing the matrix elements with the propagators 
in terms in hypergeometric functions, following Refs. [15,16]. Final values are indicated in Table [TJ 

Table 1: Numerical results for the decay rates nS -^ IS and nD -^ IS for hydrogen. 
The rates scale with Z^ for a hydrogenlike ions with nuclear charge number Z. Units 
are inverse seconds. To obtain the decay rate in Hertz, one needs to divide by a factor 
of 27r. We here supplement the results given in Ref. [1] by some values for higher 
excited S states and we also indicated results for nD ^ IS, which were not treated 
in Ref. [1]. 
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The one-loop as well as the two-loop self-energy shifts of hydrogenic states are well known to follow scaling laws of 
the form of inverse powers of the principal quantum number n, as analyzed in Ref. [17]. The two-photon decay rate is 
the imaginary part of this energy shift and is thus expected to follow an analogous trend with the principal quantum 
numbers. Analyzing the data in Table [TJ we find that the nD ^15* state results appear to follow the asymptotic 
behaviour (expressed in inverse seconds) 
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whereas for nS — > 15 decay, a fractional power apparently leads to a more satisfactory representation of the data, 
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The results indicated in Table [T] are consistent with a decrease of the two-photon decay rate with increasing n. 

4 Discussion and Comparison 

When comparing to the existing literature, it is useful, first of all, to note the calculations [2,5,6], which are apparently 
based on second-order perturbation theory for the two-photon transition amplitude. As a consequence, they present 
singularities when the energy of one of the photons reaches a level situated between the initial and final states, and no 
procedure is given in the cited references if one does not go beyond second order. When evaluating differential transition 
rates (Refs. [5, 18, 19]), the absence of the infinitesimal imaginary part does not matter, and the numerical results in 
the velocity gauge [5, 19] and in the length gauge [18] fully agree. The problem arises when one tries to evaluate the 
total decay rate, as the existing singularities are not integrable. Although in Ref. [7] fourth-order perturbation theory 
was used, a consistent answer does not appear to have been found. 

It appears that in general, two approaches have been used so far in the literature in order to deal with the problematic 
double poles for the photon energy integrations: (i) the explicit removal of particular states from the sum over virtual 
states, and (ii) the inclusion of a width for the intermediate, virtual states. 

Let us begin the discussion with the removal of states. Indeed, Chluba and Sunyaev [2], Florescu et al. [6] as well as 
Cresser et al. [7] have used different formulas than those used here, in order to evaluate the two-photon decay rates. 
In particular, they use instead of Eq. (|14p the following formula for nS -^ IS decays, 
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where N — n (Chluba and Sunyaev, Ref. [2]) or or A^ = n+ 1 (Florescu et al. [5,6] and Cresser et al., Ref. [7]), and the 
notation ly > N oi course means that one should sum over the discrete spectrum for all virtual states with principal 
quantum numbers as indicated, and of course integrate over the entire continuum spectrum in addition. For nD -^ IS 
decays, the cited authors use 
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with the same proposed values for N . In this case, because the problematic virtual states of lower energy than the initial 
state |(/>i) have been explicitly removed from the sum over virtual states, there are no more singularities infinitesimally 
displaced from the integration contours present, and there is therefore no need for any infinitesimal imaginary part ie 
in the propagator denominators. Furthermore, | • [^ is equivalent to (•)^ provided our assumption formulated in Eq. ([5]) 
holds. The corresponding velocity-gauge expressions. 
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are not equivalent to the length-gauge expressions in Eqs. (fT5)) and (PO)) . because the relation P^ breaks down if 
the sum over v does not extend over the entire hydrogen spectrum. The explicit removal of the "problematic" virtual 
states from the propagators avoids the necessity of indicating the infinitesimal imaginary terms in the propagator 
denominators, but the removal operation leads to different expressions in the length and the velocity gauges and is 
thus not gauge invariant. 



To illustrate this finding by a numerical example, we observe that we can reproduce the value of 73^ = 0.131 813 s ^ 
for the decay 3-D — > IS* with N = n+l using the length-gauge expression ((20|) . in agreement with Eq. (20) of Ref. [7]. 
However, the velocity gauge expression ([^^ gives a different result, namely rj^jj — 0.439 368 s""'^. These two results 
have to be contrasted with the gauge-invariant result of Fg^ = 1.042 896 s^^, indicated in Table [TJ For the decay 
3S -^ IS, the values are 73^ = 8.225 796 s~^ in agreement with Eq. (19) of Ref. [7], and the velocity-gauge result with 
2P and 3P virtual states removed is 773^ — 6.192 881 s^^, whereas the gauge-invariant result with the full hydrogenic 
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spectrum of virtual states reads Fgg, = 2.082 853 s ^ (see Table [T]). It is interesting to observe that Fg^ > 73^), but 
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Let us now turn our attention to the inclusion of a decay width for the intermediate states. Indeed, the authors of 
Refs. [2,5-7] arrive at the expressions ([21]) and (|22|) after analyzing the expression (for illustrative purposes we restrict 
ourselves here to the nS -^ IS decay) 
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Let us consider 35 state as an example. The only "problematic" virtual state is the 2P state [v 
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it is possible to show, that the term with v = 2 va the expression ([23]) gives rise to a contribution which is equivalent 
to the one-photon decay rate 35' -^ 2P, and this decay rate is just the total one-photon decay rate of the 35 state, 
and it is equal to the imaginary part of the one- loop self-energy of the 35 state (in the dipole approximation) . The 
authors of Refs. [2, 5-7] thus conclude that this term should be interpreted as the one-photon decay rate of the 35 
state, which has got nothing to do with the two-photon decay process, and this observation appears to be the basis 
for their removal of the 2P state from the sum over virtual states (see also the analysis in footnote 4 of Ref. [7]). 

Despite the appealing aspects of the removal operation, it is unfortunately not gauge invariant, as shown above, and 
we would like to point out two more aspects that merit a discussion. First and foremost, the discussion in footnote 4 of 
Ref. [7] shows that the expression ([23]) gives rise to a one-photon decay rate, effectively mixing the two-loop self-energy 
with the one-photon self-energy (according to the interpretation of a decay rate as an imaginary part of an energy 
shift). If one would take the expression ([23]) literally, then one should be careful to avoid a double counting of the 
one-photon decay rate, which is already contained in the one-loop self-energy and should not be obtained once more 
from the imaginary part of the two-loop self-energy. Cascade contributions are discussed in more detail in Sec.[5]below. 

The second aspect is observed when the analysis in footnote 4 of Ref. [7] is generalized to the 45 ^ 15 decay. In that 
case, two cascades are possible, namely 45 -^ 3P -^ 15 and 45 -^ 2P -^ 15. As an easy generalization of the analysis 
in footnote 4 of Ref. [7] shows, the full one-photon decay rate of the 45 state is obtained from the expression (P5)) 
only if the virtual 2P and 3P are endowed with their partial decay rates to the 15 ground state, i.e. the 3P decay 
rate should be inserted into the propagator denominators as the partial decay rate 3P -^ 15, excluding the decay 
process 3P —^ 25. If one generalizes these considerations further, namely to a general decay nS — > 15, then this would 



to regularize the divergence in l/FJ, in Eq 



imply that one should use different decay rates Fi, in Eq 
adjusting them according to the decay process under study. That prescription would be highly counterintuitive as the 
virtual states should somehow "know" about properties of the initial and final states of the decay process. The ensuing 
questions have already been noticed by Chluba and Sunyaev [2]. 

Let us conclude this section with a remark on the asymptotics pT|) and (|18p . which permit an extrapolation of our 
results to Rydberg states with high principal quantum numbers. Some investigations, including Ref. [2], lead to results 
for the two-photon decay rates of higher excited state which exhibit a linear increase with n instead of a decrease with 
at least n~^, as indicated in Eqs. p?)) and p8|) . It is well-known that the one-photon rates decrease approximately 
with n~^ (see Ref. [14]). If the two-photon rates would indeed increase linearly, then there would be a relative factor 
n^ with which two-photon rates would grow in comparison to one-photon rates as the principal quantum number 
of a state increases. If we take into account the relative scaling factor of Z'^o? j-n by which two-photon rates are 
suppressed with respect to one-photon rates, then we would have to conclude that the two-photon rates overtake the 
one-photon rates for states with a comparatively low principal quantum number of n sa 50/-\AZ in a hydrogenlike 
ion with nuclear charge number Z. For our results as indicated in Table [U the two-photon rates are suppressed with 
respect to one-photon rates by a relative factor Z'^o? ji^ for all hydrogenic states, because the scaling with n is obtained 
to be approximately the same for the one- as well as the two-photon rates, and the natural hierarchy of the likelihood 
of one- and two- photon events is preserved for all states. 



5 Extraction of the Cascade Contribution 

As in Sec. [H let us focus on a particular example whose generalization is obvious, namely (this time) the 3S -^ IS 
decays, for which the cascade 35 -^ 2P -^ IS needs to be addressed. Let us go back once more to Eq. (|14l) . 
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and adopt the cumbersome, but absolutely unique notation P.V. for the principal value part of the distribution. If we 
use the formula 
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for all propagator denominators in Eq. (|25p and extract only the contribution due to the delta functions, then the only 
contributing virtual state is the 2P state. Because there is a product of two terms both of which become singular, we 
cannot avoid to obtain the square of the delta function, 

S^{lu - E3S - E2P) = 5(0) 5{uj - E3S - E2P) = ^5{oj- E^s - E2p) , (27) 

zvr 

and a further term proportional to S{lu — E2P + -Eis) T/(2tt). Here, T is the (long) observation time proportional to 
(5(0) in energy space (see, e.g. Ref. [3]). The sum of the terms proportional to (5(0) reads 

'-^3S ~ ^ ^ 3S^2P^ 2P^1S ' / ^^ ^3S ~ n ^ 3S^2P ^ 2P^15 ' \'^°) 

where we introduce an obvious notation for the partial one-photon rates ^\,g^')p ^^^d ^2P^is- Note, in particular, 
that the resulting expression for Cgg is gauge invariant. 

Our result (P5|) has just the right form to describe the cascade decay, except for the "wrong" sign. For the term to 
contribute to the decay of the ZS state, it should be positive, but it turns out to be negative. Let us defer a discussion 
of this issue and instead consider the extraction of the cascade contribution from the expression 

f^ = ^ hm ^TL .3 (^3, - ^,, - .)3 TT V ( (1^ 1 1-1 1 ^P) (-^ 1 1-1 1 3^) ^ (1^ I \A I ^P) i^P 1 1x1 1 35) | 
3S 27^m2 e^o J ^ 3S IS ^ 1 i Z^ \ E^s ~E,p-w±it Eis -E,p+uj±ie J ' 

± I' 

(29) 
where Y[± nieans the product of two terms, with either sign. The product over the two terms with ±ie is of course 
equivalent to the square of the modulus of the two terms in the integrand, in analogy to Eq. ((23|) . If we now use ([26]), 
then we obtain 

^(2)_^^(i) p(i) /"HTr^(2) _ ^p(i) p(i) /on^ 

^3S — ^ ^ 3S-.2P ^ 2P-+1S ' / ^^^ "-'3S ^ 2 3S^2P ^ 2P— IS ■ V'J^J 

This result has the "right" sign, and it has, for large T, the right temporal dependence for a cascade process. 

In order to resolve the paradox, one first should notice that both signs found in Eqs. (j28|l and ([30]) actually have a valid 
interpretation. The two-loop self-energy contains both radiative corrections to the one-photon decay as well as the full 
two-photon decay amplitude. The radiative corrections to one-photon decay are obtained by "cutting" appropriate 
internal lines in the diagrams, and indeed, we can rederive the first three radiative corrections to one-photon as given 
in Eq. (27) of Ref. [20] by considering resonant intermediate states in the "outer" electron propagators in the terms of 
Eq. ([3]). (The remaining terms used in Eq. (27) of Ref. [20] follow from standard third-order perturbation theory.) The 
magnitude of the radiative corrections to one-photon decay is decreased by the possibility of cascade processes, due to 
the virtual-to-real conversion of the photons appearing in the integrals for the radiative corrections to the one photon 
decay at resonance, and this decrease is consistent with the sign of the right-hand side of Eq. ((28l) . On the other hand, 
the two-photon decay amplitude should be increased by the cascade processes, and this increase is consistent with the 
sign of the right-hand side of Eq. ([30]) . 

In the sum of the radiative corrections to one-photon decay and the two-photon decay, the incoherent cascade contri- 
butions cancel, and this is in analogy to the discussion in Ref. [3] for a different, but physically related process, namely 
the coherent/incoherent pair production via a virtual/real photon intermediate state by an electron in crossed, static 
electromagnetic fields. 

Immediately, new questions arise. Our considerations suggest that our formulation in Eq. Q provides infinitesimal 
imaginary parts that are appropriate for the evaluation of radiative corrections to the one-photon decay, but provides 



the "wrong" sign of the cascades for two-photon decays. This could lead to new doubt regarding whether we can extract 
a valid expression for the two-photon decay rate from our Eq. ([3]) in the first place. The question is: Can we extract, 
by some mathematically justifiable procedure, from Eq. ()29p . an expression for the two-photon decay rate which either 
confirms or invalidates our result for the two-photon decay rate, under a suitable gauge-invariant subtraction of the 
cascade contribution from the integrand in (j29p ? 

First, since the cascade contributions correspond to the delta function in Eq. ([26| . it is clear that the two-photon 
decay rate corresponds to the product of two principal- value distributions of the form. 
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(31) 



which is integrated over u. As similar problems have occurred in field theory [see Eq. (6.23) on p. 168 of Ref. [4]], we 
are provided with a guiding principle for the calculation. Namely, we consider an arbitrary function /, integrated over 
a finite interval (0,Ci;max) with /(O) = /(wmax) = 0: 
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This subtraction, applied to Eq. (HH), gives rise to 
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(32) 
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Here, we have subtract the cascade-generating terms according to Eq. ([32)l . thus leading to an integral which is finite 
under a principal- value prescription, because the double poles have explicitly been subtracted. Because the prescription 
P2p takes the numerators to exact resonance, the subtraction terms in p3p are gauge invariant and indeed proportional 
to r35^2P ^ip^is- "^^^ additional term T is due to the boundary term found in (|32p . 



T^- 2{Ess - Eis) {Ess ~ E^p)^ {E^p - E^s)^ (15 ||f || 2P)' (2P ||f|| 35)' 
Finally, returning to our original ie prescription, we have according to Eq. ()16p . 
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and in view of ([32]) and (|35p . we obtain the (perhaps somewhat surprising) equality 
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which is subject to the interpretation of the squared principal- value contribution according to Eq. (|32p . We can finally 
state that the result ([33]) agrees with formula (|25|) . so that, under the provisions of the regularization implied by 
Eq. ()32p . it is irrelevant if we start from an expression where the integrand for the two-photon decay is formulated as 
a modulus squared or with two infinitesimal imaginary parts "pointing in the same direction." 



6 Conclusions 

We have analyzed two-photon decay processes involving nS — > IS' and nD -^ IS channels in hydrogenlike ions. Our 
general formulas (j4|) and (jlip are gauge-invariant and are obtained with otherwise unspecified, arbitrary infinitesimal 
imaginary parts i e, provided the limit e — > is taken after the integrations over the photon energies have been 
performed (non-uniform convergence). Numerical results are presented in Table [1] These are nonrelativistic results 
which scale as Z^ with the nuclear charge number Z. For a relativistic generalization, see [21]. 

From a more philosophical point of view, we can say that the two-photon decay process turns out to be an extremely 
subtle physical phenomenon, which demands a lot of mathematical sophistication in its analysis. Without a careful 
handling of the distributions, including ill-defined squares of delta functions, it is impossible to obtain consistent 
answers. The current work attempts to provide a proposal for a consistent framework in which the resonant intermediate 
states and the generated double poles can be addressed, while preserving the interpretation of the integrand of the 
two-photon decay rate as a differential decay rate with respect to the photon energy. 

Three remarks conclude this work, (i) Following Ref. [3], we should point out that there is no guarantee that the 
coherent two-photon decay rate as evaluated here always needs to be positive (except for the 25* state, where no 
resonant intermediate states are present). Indeed, as Eq. ([55t shows, the result is obtained as a subtracted integral, 
and the integrand is not necessarily positive. For all transitions considered here, the rate is positive (see Table [T]), but 
it is known that radiative corrections to decay rates can be negative, and the coherent two-particle contribution to 
a decay rate beyond the cascade constitutes a correction to a decay rate which need not be positive. This statement 
is paradoxical, but we can point out that this statement has already been confirmed after Eq. (20) of Ref. [3] in an 
absolutely analogous situation, (ii) The observation time T as implied in Eq. (PD|) has to be sufficiently large (larger 
than the typical formation time of radiation in the system, according to Ref. [3], or otherwise the decay process will 
proceed in a different way). In our case, the natural formation time of radiation is given by a time inversely related to 
the decay width of the decaying states, which is naturally identified as the one-photon decay rate of the highly excited 
states, (iii) It may seem that the agreement of the integration around the infinitesimally displaced poles as described in 
Sec. [Hand the regularized principal- value prescription as described in Sec. [5] is purely accidental. However, one should 
remember that similar integrals appear in Lamb-shift related self-energy calculations (e.g., Ref. [22]), and therefore, 
the predictions for the Lamb shift of excited states would have had to be reinvestigated if we had not found agreement 
of the two computational schemes discussed here. Fortunately, the internal consistency of mathematics protects us 
from having to reinvestigate accurate theoretical predictions based on quantum electrodynamics. 

Acknowledgment 

Helpful discussion with A. Surzhykov, A. I. Milstein, P. J. Mohr, F. Khan and K. Pachucki are gratefully acknowledged. 
This work was supported by the Deutsche Forschungsgemeinschaft (Heisenberg program, contract JE285/3-1). 

References 

[1] U. D. Jentschura, J. Phys. A 40, F223 (2007). 

[2] J. Chluba and R. A. Sunyaev, e-print 0705.3033 [astro-ph]. 

[3] V. I. Ritus, Nucl. Phys. B 44, 236 (1972). 

[4] V. N. Baier, V. M. Katkov, and V. M. Strakhovenko, Electromagnetic Processes At High Energies in Oriented 
Single Crystals (World Scientific, Singapore, 1998). 

[5] V. Florescu, Phys. Rev. A 30, 2441 (1984). 

[6] V. Florescu, I. Schneider, and I. N. Mihailescu, Phys. Rev. A 38, 2189 (1988). 

[7] J. D. Cresser, A. Z. Tang, G. J. Salamo, and F. T. Chan, Phys. Rev. A 33, 1677 (1986). 

[8] K. Pachucki, Phys. Rev. A 63, 042503 (2001). 

[9] K. Pachucki and U. D. Jentschura, Phys. Rev. Lett. 91, 113005 (2003). 

[10] R. Barbieri and J. Sucher, Nucl. Phys. B 134, 155 (1978). 



[11] L. D. Landau and E. M. Lifshitz, Quantum Mechanics (Volume 3 of the Course of Theoretical Physics) (Pcrgamon 
Press, London, 1958). 

[12] D. A. Varshalovich, A. N. Moskalev, and V. K. Khcrsonskii, Quantum Theory of Angular Momentum (World 
Scientific, Singapore, 1988). 

[13] M. Goppert-Mayer, Ann. Phys. (Leipzig) 9, 273 (1931). 

[14] H. A. Bcthe and E. E. Salpeter, Quantum Mechanics of One- and Two-Electron Atoms (Springer, Berlin, 1957). 

[15] M. Gavrila and A. Costcscu, Phys. Rev. A 2, 1752 (1970). 

[16] K. Pachucki, Ann. Phys. (N.Y.) 226, 1 (1993). 

[17] U. D. Jcntschura, E.-O. Le Bigot, P. J. Mohr, P. Indehcato, and G. Soff, Phys. Rev. Lett. 90, 163001 (2003). 

[18] J. H. Tung, X. M. Ye, G. J. Salamo, and F. T. Chan, Phys. Rev. A 30, 1175 (1984). 

[19] V. Florescu, S. Patrascu, and O. Stoican, Phys. Rev. A 36, 2155 (1987). 

[20] J. Sapirstein, K. Pachucki, and K. T. Cheng, Phys. Rev. A 69, 022113 (2004). 

[21] A. Surzhykov and U. D. Jcntschura, submitted (2008). 

[22] U. D. Jcntschura, G. Soff, and P. J. Mohr, Phys. Rev. A 56, 1739 (1997). 



10 



